Spin Nernst effect and Nernst effect in two-dimensional electron systems 
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We study the Nernst effect and the spin Nernst effect, that a longitudinal thermal gradient 
induces a transverse voltage and a spin current. A mesoscopic four-terminal cross-bar device having 
the Rashba spin-orbit interaction (SOI) under a perpendicular magnetic field is considered. For zero 
SOI, the Nernst coefficient peaks when the Fermi level crosses the Landau Levels. In the presence of 
the SOf, the Nernst peaks split, and the spin Nernst effect appears and exhibits a series of oscillatory 
structures. The larger SOI is or the weaker magnetic field is, the more pronounced the spin Nernst 
effect is. The results also show that the Nernst and spin Nernst coefficients are sensitive to the 
detailed characteristics of the sample and the contacts. In addition, the Nernst effect is found to 
survive in strong disorder than the spin Nernst effect does. 

PACS numbers: 72.15.Jf, 72.25.-b, 73.23.-b, 73.43.-f 



The Hall-like effect, namely, a longitudinal force in- 
duces a transverse current, has been a fascinating topic 
since the early days of the condensed matter physics. The 
integer and the fractional quantized Hall effects, two cel- 
ebrating examples, have been extensively investigated for 
the last three decades, but remain to be active research 
fields till now. Recently, another Hall-like effect, spin 
Hall effect, in which a longitudinal voltage bias induces a 
transverse spin current due to the scattering by magnetic 
impurities or due to the existence of a spin-orbit inter- 
action (SOI), has generated a great deal of interest.--- 
Apart from a large number of theoretical studies, several 
experimental investigations also made important contri- 
butions to the field. Up to now, SOI has indeed been 
found to be substantial in some semiconductors despite 
it is a relativistic effect, and its strength can be tuned by 
the gate voltage in the experiment i^i^ In particular, the 
spin Hall effect has been detected experimentally by ob- 
serving the transverse opposite-spin accumulations near 
the two edges of the sample i^i^ 

The Nernst effect, a thermoelectric property, in which 
a longitudinal thermal gradient induces a transverse cur- 
rent (or a bias ISV with open boundary) while under 
a perpendicular magnetic field, is also a Hall-like effect. 
The thermoelectric coefficients (including the Seebeck co- 
efficient and the Nernst coefficient) of electronic systems 
are known to be more sensitive to the details of the den- 
sity of states than the conductancCjMA and these de- 
tailed information of the density of states is importance 
for the design of the electronic devices. But the ther- 
moelectric measurement is usually more difficult to carry 
out than the conventional transport measurements, par- 
ticularly for low-dimensional systems or nano-devices. 
Fortunately, because of the development of the micro- 
fabrication technology and the low-temperature measure- 
ment technology in the last two decades, the thermoelec- 
tric measurement in low-dimensional samples has been 
feasible nowj ^^d^ Recently, the thermopower of the quan- 
tum dot was measured, and the results in the Kondo 
regime show a clear deviation from the semiclassical Mott 
relatioufiSii^ The Nernst effect in bismuth has also been 



detected and the Nernst coefficient peaks at positions 
when Fermi level crosses over the Landau levels (LLs) 
Meanwhile, limited theoretical studies of the Nernst ef- 
fect have also appearedfi^' 

In this paper, we study the Nernst effect and spin 
Nernst effect in a two-dimensional electron gas with a 
SOI and under a perpendicular magnetic field B. For 
the first time, the spin Nernst effect, a novel Hall-like ef- 
fect, is investigated. The spin Nernest effect implies that 
a longitudinal thermal gradient AT induces a transverse 
spin current. The spin Nernst coefficient should be more 
sensitive to the details of the spin density of states of the 
system than the spin Hall conductance, similar as their 
electronic counterparts. ^'^'^ We consider the system as 
shown in Fig. la, consisting of a square center region con- 
nected to four ideal semi-infinite leads. A longitudinal 
thermal gradient AT is added between the leads 1 and 3. 
This thermal gradient induces a transverse Hall voltage 
Vh with the open boundary condition under a perpen- 
dicular magnetic field _B; a transverse spin current Jsh 
in the closed boundary condition with a SOI. By using 
a tight-binding model and the Landauar-Buttiker (LB) 
formula with the aid of the Green's functions, the Nernst 
coefficient {N^ = Vh/^T) and spin Nernst coeffi- 
cient Ns {Ns = Jsh/ ^T) are calculated. Without a SOI, 
the Nernst coefficient peaks when the Fermi level Ep 
crosses the LLs, and spin Nernst coefficient Ng is absent, 
consistent with the recent experimental findings<i2. In the 
presence of a SOI, each LL splits into two, consequently, 
each Nernst peak splits into two peaks. Meanwhile, the 
spin Nernst effect emerges and its coefficient Ng exhibits 
a series of oscillatory structures. The oscillation is en- 
hanced with increasing SOI but is damped by a large B. 
In addition, the Nernst effect is found to survive in strong 
disorder than the spin Nernst effect does. 

In the tight-binding representation, the cross-bar sam- 
ple is described by the Hamiltonian>i^, 
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FIG. 1: (color online) (a) Schematic diagram for the four- 
terminal cross-bar sample with a thermal gradient AT applied 
between the longitudinal lead-1 and lead-3. (b) Schematic 
view of the LLs of the center region, the Fermi energy Ef, and 
the bias V2. The oscillatory line across Ep is the difference of 
ME)~f2iE). 



- cls.J^'^yU'^i<r'^''"''+H.c.] (1) 

where cl^{cicr) is the creation (annihilation) operator of 
electrons in the site i = (n, m) with spin a. t — 
fi^ /{2m*a?) is the hopping matrix element with the lat- 
tice constant a, and 5x and 5y are the unit vectors along 
the X and y directions. £i is the on-site energy, which 
is set to everywhere for the clean system. While in a 
disorder system, e\ in the center region is set by a uni- 
form random distribution [-W/2,W/2]. The last term 
in Eq.(l) represents the Rashba SOI with being its 
strength. In order to avoid confusion in calculating the 
spin current, Vr is set to zero in the lead- 2 and lead-4. 
The extra phase 6 — ea?B/h is from the perpendicular 
magnetic field B. Here the Zeeman effect and electron- 
electron interaction are neglectedf^^ The Zeeman split 
could be small in some of the two-dimensional electron 
systems. The electron-electron interaction is weak in sys- 
tems with high carrier density. 

The particle current Jpu in the transverse lead-p with 
spin (7 =1: J. can be obtained by the LB formula>i^ 



J DfT 
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dE Tp,^,{E)[f,{E) - f,{E)] (2) 



where Tp„ q{E) is the transmission coefficient from the 
lead-q to the lead-p with spin a and E is the energy of 
the incident electron. The transmission coefhcient can 
be calculated from: Tp^,q{E) = TrlTp^G'TgG"], where 
the line-width function Tp,, (i;) = i(S;;^ - S;;t), Tg = 
+ Tqi, and S^^ is the retarded self-energy due to 
coupling to the lead-p with spin a. The Green's function 
G-iE) = [G-{E)]^ = {EI-Uo-Ep. and Hq is 

the Hamiltonian of the central region. fp{E) in Eq.(2) is 
the electronic Fermi distribution function of the lead-p, 
and fp{E) = l/{exp[{E - Ep - Vp)/kBTp] + 1} with the 
bias Vp and temperature Tp. After getting the particle 
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FIG. 2: (color online) Ne (black) and Na (red or gray) vs. 
1/B for different Vr = (a), 0.02 (b), 0.05 (c), and 0.1 (d). 
The other parameters are Ep ~ —3, the size of the center 
region L = 40a, and T — 0.01. In the inset of (a), the dots 
are the inverse of peaks' maxima vs. TVth peak, and the line 
is {N+^)/ln2 vs. iV. 



current Jp^, the (charge) current is Jj 



pe 



+ Jpl) 



e(Jp^ 

and the spin current is Jps = (?i/2)(Jp| — Jpi). 

Considering a small temperature gradient AT and zero 
bias applied on the longitudinal lead-1, 3, we can set the 
temperatures Ti = T +AT /2, T3 = T-AT /2, T2 = Tt 
T, and the biases Vi = V3 = 0. From the open boundary 
condition with J2e = J-ie = 0, the transverse voltage V2 
and V4 can be obtained, and consequently the Nernst 
coefficient A^e = (V2 — V4)/AT. In the clean system, Ne 
is expressed as 



1 jdE {T2i~T23){E~Ep)f{l-f) 
sT / di; (Tai + T23 + 2T24)/(1 - /) ' 



(3) 



where — T2]^p + ?2X,p- The spin Hall current J2S and 
J4S are calculated with the closed boundary condition 
having V2 — V4, = 0. In the clean system, J2s — — J is be- 
cause of the symmetry property of the system. In a dirty 
system, J2S may not equal to — J is for a given disorder 
configuration, but J2S — —J is still holds after average 
over many configurations. The spin Nernst coefficient 
Ns = J2s/AT, and can be reduced to: 

where AT2p = ?2T,p ^ ^2j.,p- At low temperature limit 
(T 0), the Nernst coefficient Ne and the spin Nernst 
coefficient usually depend linearly on temperature. 

But while Tutp^^7f-t\e=Ep = ( ^(ATai-ATaa) b^-Ef = 

- T23 {AT21 - AT23) at E = Ef 



d(T2i-T2: , 

0), or in other words T21 - 
is discontinuous, A^e (-^s) is temperature independent. 

In the numerical calculations, we set t = ?iV(2w*a^) 
as the energy unit and ^a^ as the unit of the mag- 
netic field B. If taking the effective electron mass 
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FIG. 3: (color online) (up panel) The LLs En vs. 1/B for 
Vr — (black) and 0.1 (red or gray). The blue dotted line 
is the Fermi level Ep- (down panel) Ne vs. 1/B for the 
Rashba SOI in lead-2, 4 is (black) and 0.1 (red or gray). The 
parameters are Ep ~ 1.0, Vr = 0.1, L = 20a, and T — 0.01. 
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FIG. 4: (color online) (black) and Ns (red or gray) vs. 1/B 
for different Fermi level Ep with the parameters Vr — 0.05, 
T = 0.01, and L = 40a. 



m* = 0.05me and the lattice constant a — 12.5nm, t 
is about 5meV, B — 1 corresponds to 4.2 Tesla, and 
Vr = Q.lt corresponds to the Rashba SOI parameter 
a = 1.25 X 10~^^cV-ni which can be experimentally mod- 
ulated by the gate voltage. We consider square samples 
and the center-region size is either L = 40a or L = 20a 
in our calculations. Temperature is fixed at T = 0.01, 
that is about IK. 

Fig. 2 shows the Nernst coefficient and spin Nernst 
coefficient A^s versus the inverse of magnetic field 1/B 
for the different SOI strength Vu in the clean system 
(W = 0). While without the SOI {Vr = 0),Ns is exactly 
zero, but Ne exhibits a series of equal spacing peaks. Ne 
peaks when the Fermi level Ep crosses over LLs, and 
it is damped when Ep lies between adjacent LLs. The 
peak interval is eh/{Tn*Ep) where Ep — Ep + 4i is the 
distance from the Fermi energy to the band bottom —At. 
The inverse of the height He of the iVth peak is linearly 
dependent on N, with He (x N + 1/2 (as shown in inset 
of Fig. 2a). Let us explain these characteristics with aids 
of the physical picture in Fig. lb. Under a strong B, the 
transmission coefficients T27i{E) and T2i{E) are usually 
zero, and T2i{E) is an integer. Then the current J2e — 
(e/h) J dET2i{E)[f2{E)- h{E)] fromEq.(2). Due to the 
thermal gradient, /2 — /i exhibits an oscillatory structure 
around Ep as shown in Fig. lb, and the electrons with 
energy above and below Ep contribute opposite signs to 
the thermocurrent J2e- When all LLs are far from Ep, 
T21 is a constant near Ep, the currents flowing in or out 
cancel each other, leading to J2e = at V2 = 0. On the 
other hand, when N LLs are below Ep but one LL is at 
Ep (see Fig. lb), a net current J2e is induced at V2 = 0. In 
the open circuit case, V2 has to be raised to make J2e = 0, 
and V2 is the ratio 1/{2N + 1). In fact, from Eq.(2) and 



assuming that T32 = T42 — and T21 is an integer for 
large B, we can analytically obtain that the peak height 
of the Nernst coefficient is He = ^In2/{N + 1/2). This 
result is identical with the result of the thermopower in 
a two-terminal systemj^^ 

While in the presence of a SOI (Vr 0), the LLs 
split. As a result, the peaks of the Nernst coefficient Ne 
also split and the spin Nernst coefficient Ng emerges (see 
Fig. 2). The splitting is more pronounced for stronger 
SOI Vji or weaker magnetic field B. The positions of the 
right sub-peaks of Ne are consistent with LLs, but not 
the left sub-peaks. To see this, we magnify the second 
peak of Fig. 2d, and also plot in Fig. 3 the LLs versus 1/B 
without SOI (Vr = 0) and with SOI {Vr ^0). It clearly 
shows that the left sub-peak is in line with the original 
un-split LL at Vr = (see mark b in Fig. 3), not in align- 
ment with the split LLs. In order to thoroughly study the 
peak positions, we also plot Ne for the uniform system, 
in which SOI exists in all parts, including the leads-2, 4. 
Now the left sub-peak moves to align with the split LL. 
(see mark a in Fig. 3). So the counterintuitive phenomena 
entirely comes from the non-uniformity of SOI, in which 
the SOI is absent in the leads-2, 4 and an interface be- 
tween Vr — and Vr ^ emerges. This interface causes 
additional scattering for an incident electron, and one of 
the edge states goes directly from lead-1 to lead-3 instead 
of from lead-1 to lead-2, so the left sub-peak position in 
Ne is moved. This means that the Nernst effect can re- 
flect the detailed structure of the transverse leads and its 
contact to the sample. This is essential difference to the 
regular Hall eSect^'^ 

Next, we study the spin Nernst coefficient Ng, which 
emerges with Vr (see Fig.2b-d). In the vicinity of the 
right sub-peak of Ne, Ng exhibits an oscillatory structure 
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FIG. 5: iVe (a) and A^^ (b) vs. 1/B for the different disorder 
strengths W = ^ (solid curve), 0.1 (dotted curve), 0.5 (dashed 
curve), and 1.0 (dash-dotted curve). The other parameters 
are Ef = -3, Vr = 0.1, T = 0.01, and L = 20a. 

and N, oc d[lnNe{EF)]l dEp. This relation of N, and 
is similar to the semiclassical Mott relation between the 
thermopower and conductanceJ^ However, Ns is quite 
small and docs not show an oscillatory structure around 
the left sub-peak of , and the Mott-like relation breaks 
there. We can qualitatively analyze these phenomena 
using Schrodinger equation, in which we can analytically 
obtain the split LLs and corresponding wave functions. It 
is found that for weak SOI (such as Vr < O.lt), the wave 
functions of the high sub-LLs are strongly spin polarized 
in the z-direction, while one of the low sub-LLs are hardly 
spin polarized. As a result, the spin current J2,4:s with 
an oscillatory structure only exists when Ep crosses over 
the high sub-LLs, which corresponds to the positions of 
the right sub-peaks of Ng. In addition, a larger SOI or a 
weaker magnetic field will cause a stronger spin Nernst 
signals (see Fig. 2), due to the competition between the 
magnetic field and the SOI. In fact, the Rashba SOI is to 
drive the electrons with opposite spins to opposite direc- 
tions transversely which leads to the spin Nernst effect, 
but a magnetic field B is to drive all electrons in the same 
transverse direction. Thus, B weakens the spin Nernst 
effect. 

The relation of the LLs E^ with the magnetic field B 
is: En ^ ^{N + l/2)~ so the period in the Nernst 
signal is e'h/{m*Ep), changeable by adjusting the Fermi 
level Ep. In Fig. 4, we plot A^e and Ng as functions of the 
inverse of B for different Ep. The results for different Ep 
show similar behaviors. With increasing of Ep, the peaks 
of Ne and the oscillatory structures of Ns are getting 



closer, and the magnitudes of iVg and Ns signals at same 
filling factors are getting weaker. 

Finally we discuss the disorder effect on the Nernst 
and spin Nernst effects. Fig. 5 shows Ne and Ns versus 
the inverse B for different disorder strength W. Here Ne 
and Ns are averaged over 500 disorder configurations. 
For a small disorder (e.g. W = 0.1), both Ne and Ns 
are hardly affected. For an intermediate disorder, such 
as W = 0.5, the two sub-peak heights of Ne almost keep 
their strengths as for W = 0, but the valley between two 
sub-peaks is greatly deepened, so that the two sub-peak 
structure is even clearer (see Fig. 5a). With further in- 
creasing of the disorder W, the left sub-peak of Ne is 
decreased while the right sub-peak is less affected, mean- 
while the oscillatory structure of Ng is weakened. How- 
ever, in the vicinity of the left sub-peak of iVe, the spin 
Nernst coefficient Ng, which is very small at = 0, is 
enhanced by W (see Fig. 5b). Finally, for very large dis- 
order W , the system goes into an insulating regime, both 
Ne and Ns vanish. 

Before summary, we would like to make a couple 
comments concerning the novel spin Nernst effect, (i) 
The spin Nernst effect is NOT a simple combination 
of the Seebeck effect and the spin Hall effect. In 
fact, the Seebeck coefficient is mainly determined by 
{dTi3{E)) / {dE)\p=Ep , while the spin Nernst coefficient 
depends on AT23 and AT21. (ii) The spin Nernst effect 
can be measured in similar ways that the spin Hall effect 
is observed;^ e.g. through spin accumulations. 

In summary, the Nernst effect and spin Nernst effect 
in a two-dimensional cross-bar with a spin-orbit inter- 
action and under a perpendicular magnetic field are in- 
vestigated. The Nernst signal exhibits a series of peaks, 
and the inverse of a peak height goes linearly to the se- 
quence number of the peak. While in the presence of a 
SOI, these Nernst peaks split, and the spin Nernst effect 
appears, which exhibits an oscillatory structure versus 
the magnetic field. The relation of the Nernst and spin 
Nernst coefficients is similar to the semiclassical Mott 
relation around one sub-peak, but has a great discrep- 
ancy around the other sub-peak. In addition, the disor- 
der effect on the Nernst and spin Nernst effects is also 
discussed. 
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